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统 的 优化 算法 往往 依赖 系统 的 精确 数学 模型 和 系统 参数 ， 对 言 息 ， 研 究 者 设计 了 一 些 新 的 优化 算法 。 例 如 教学 优化 算法 


基于 多 子 群 的 社会 群体 优化 算法 
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摘 要 : 社会 群体 优化 算法 (social group optimization，SGO) 是 一 种 基于 社会 群体 学 习 而 提出 的 一 种 新 型 优化 算法 。 针 
对 社会 群体 优化 算法 易于 陷入 局 部 最 优 问题 ， 提 出 了 一 种 多 群 社会 群体 学 习 算法 (MPSGO)。 本 算法 采用 多 子 群 学 习 
方法 ， 对 算法 两 个 阶段 的 个 体 学 习 方法 进行 改进 ， 在 维持 群体 收 你 性能 的 前 提 下 提高 群体 多 样 性 ， 同 时 对 部 分 个 体 中 
引入 量子 学 习 ， 使 个 体 学 习 的 有 用 信息 得 以 增强 ; 此 外 ， 每 隔 一 定 代数 对 子 群 进行 随机 重组 ， 既 能 保证 各 子 群 个 体 充 
分 进化 , 又 维持 了 子 群 多 样 性 。 在 设计 算法 的 基础 上 , 分 析 了 其 收敛 性 和 多 样 性 ; 通过 与 其 他 四 种 算法 进行 对 比 实 验 ， 
验证 了 改进 后 算法 性 能 更 优 。 
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Multiple subgroups based social group optimization algorithm 


Liu Yajun, Chen Debaoi, Zou Feng, Li Zheng, Wang Suxia 
(School of Physics & Electronic Information, Huaibei Normal University, Huaibei Anhui 235000, China) 


Abstract: Social group optimization (SGO) is a novel optimization algorithm which is based on social group learning. This 
paper proposed a multi-group social group learning algorithm (MPSGO) to solve the problem that the social group optimization 
algorithm was easy convergent to local optima. This algorithm adopted the multi-subgroup learning method, and the 
improvement of the individual learning method in the two stages of the algorithm. The diversity of the population was improved 
with maintaining the convergence of the population. At the same time, it introduced quantum-behaved learning method into 
MPSGO for parts of individuals to enhance the useful information of individual learning. In addition, the population was 
randomly regrouped to generate new subgroups after a certain generation. The diversity of subgroups wais maintained and the 
individuals in every subgroup were fully evolved. Based on the designing of algorithm, this paper analyzed the algorithm of 
convergence and diversity. Compared with the other 4 algorithms, it is proved that the performance of the improved algorithm 
is better. 
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这 些 算 法 都 利用 个 体 间 信息 的 交换 或 产生 部 分 新 信息 实现 对 问 
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优化 问题 是 现实 生活 中 经 常 遇 到 的 问题 ， 如 何 设计 算法 | 


S 


K 数 往往 难以 确定 ， 一 定 程 度 上 影响 了 算法 的 应 用 。 


速 准 确 地 得 到 问题 的 最 优 解 ， 是 解决 优化 问题 的 关键 一 环 。 传 近年 来 ， 为 减少 算法 自身 待 确定 参数 数量 ， 充 分 利用 群体 


题 的 优化 ， 算 法 自身 参数 对 算法 性 能 影响 较 大 ， 且 最 优化 的 参 


些 模 型 难以 建立 的 系统 ， 传 统 优化 方法 难以 解决 。 自 然 计算 方 。 (TLBO) 申 ， 模 拟 班级 教师 的 授课 和 学 生 学 习 过 程 ， 教 师 在 班级 


法 模拟 生物 的 自然 或 社会 特性 ， 其 不 依赖 系统 的 精确 模型 ， 已。 分 享 自己 的 知识 ,提高 班级 的 


六 


T 


均 成 绩 , 学生 除 向 教师 学 习 外 ， 


得 到 广泛 应 用 。 经 典 的 自然 计算 方法 ， 如 遗传 算法 (GA)m、 粒 还 向 班级 的 其 他 学 生 学 习 ， 以 进一步 提高 自己 的 成 绩 。 此 方法 
子 群 算法 (PSO) 中 、\ 差分 进化 算法 (DE)BJ、 进化 策略 (ES) 外 等 都 已 。 中 由 于 引入 教师 的 信息 ， 算 法 收敛 速度 较 快 ， 但 由 于 个 体 每 次 
成 功 应 用 于 函数 优化 、 工 程 设 计 、 自动 控制 ,信号 处 理 等 领域 。 学 习 都 保留 最 好 的 解 ， 群 体 多 样 性 丢失 较 快 ， 算 法 易于 陷入 局 
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部 最 优 解 。 步 提 高 TLBO 算法 的 性 能 ， 研 究 者 提出 了 一 
些 改进 算法 。 例 如 基于 邻 域 搜索 的 教学 优化 算法 (NSTLBO)I9 在 
教学 过 程 中 引入 个 体 的 邻 域 信息 ， 以 增加 个 体 学 习 方 法 的 多 样 
性 ， 提 高 了 TLBO 算法 的 性 能 。 精 英 教 学 算法 (ETLBO IE 


其 中 : c 为 自我 反省 参数 ， 其 取 值 通常 在 0-1 间 ， 其 经 验 值 外 
7402; rJ 0~1 间 的 随机 数 ， &best 为 当前 代 最 优 个 体 的 第 j 
基因 Xnew,, 与 Xold,, 分 别 为 第 i 个 个 体 更 新 前 后 的 第 j 
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因 。 


化 过 程 中 保留 多 个 精英 个 体 ， 利 用 精英 个 体 的 引导 作用 ， 改 变 
个 体 的 跟踪 行为 ， 进 一 步 提高 算法 的 收敛 精度 。 

为 进一步 提高 算法 的 收敛 速度 , 充分 发 挥 个 体 的 学 习 能 力 ， 
Satapathy 等 人 [8] 于 2016 年 提出 的 一 种 基于 社会 群体 学 习 的 新 
型 优化 算法 (SGO)。 此 算法 操作 简单 ， 易 于 实现 ,在 一 些 问题 
优化 中 取得 了 良好 的 效果 外。 此 算法 将 个 体 的 优化 过 程 分 为 两 
个 阶段 ， 即 提高 阶段 和 获得 阶段 。 在 提高 阶段 ， 个 体 以 群体 当 
前 最 好 解 作 为 学 习 向 导 ; 在 获得 阶段 ， 个 体 以 群体 中 任意 一 个 
个 体 和 当前 代 群 体 最 好 解 作为 指导 进行 学 习 ， 算 法 的 收敛 速度 
较 快 。 一 些 典 型 数据 集 上 的 测试 结果 显示 了 SGO 算法 的 有 效 


多 
维 
1.3 获得 阶段 
在 获得 阶段 ， 每 个 个 体 从 种 群 中 随机 选择 一 个 个 体 作 为 学 
习 对 象 ， 并 结合 当前 代 群 体 的 最 优 个 体 信 息 进 行 学 习 。 县 体 学 
习 方 法 如 式 (3) (4) 所 示 。 
if f(X,) isbetter than f(X.) 


Xnew, , = Xold, ;+ *(x,, -Xij )*n *(gbest, -X,) (3) 


else Xnew,, = Xold, ;+ *( X , -X,)*n *(gbest, -X,) (4) 


其 中 : X 与 Xj 分 别 为 个 体 i 和 个 体 的 第 j 维基 因 ; Xnew,; 


PEB, BAR SGO 算法 的 收敛 速度 较 快 ， 但 其 至 少 存在 两 方面 
的 不 足 : 首先 ，SGO 算法 在 学 习 的 两 个 阶段 中 均 以 群体 最 好 个 
体 作为 向 导 ， 算 法 的 收敛 速度 较 快 ， 但 群体 多 样 性 容易 丢失 ， 
导致 算法 易于 陷入 局 部 最 优 解 ; 其 次 ,由 于 算法 学 习 方法 单一 ， 
旦 陷入 局 部 最 优 解 ,个 体 逃 离 局 部 最 优 解 的 能 力 不 足 .因此 ， 
在 保持 SGO 算法 收敛 速度 的 前 提 下 提高 群体 的 多 样 性 ， 成 为 
提高 算法 整体 性 能 的 关键 。 本 文 在 对 个 体 更 新 方程 改进 的 基础 
上 ， 利 用 多 群 算法 具有 维持 群体 多 样 性 的 优点 ， 结 合 量子 学 习 
强化 所 需 解 的 能 力 ， 设 计 了 一 种 多 群 社会 优化 算法 。 本 文 的 贡 
献 主要 有 三 点 : a) 引 入 多 子 群 学 习 提 高 群体 的 多 样 性 ; b) 对 每 个 
子 群 中 的 个 体 ， 在 获得 阶段 随机 引入 其 他 子 群 最 优 解 ， 并 将 其 


与 Xold,, 分 别 为 第 i 个 个 体 更 新 前 后 的 第 j 维基 因 ; abest; 为 
当前 代 最 优 个 体 的 第 j 维基 因 ; n. n AA 0-1 间 的 随机 数 。 


TT 


2 多 子 群 SGO 算法 (MPSGO) 


在 基本 SGO 算法 的 两 个 阶段 中 ， 个 体 均 以 当前 代 群 体 的 
最 优 个 体 为 参考 进行 位 置 更 新 。 实 践 表明 ， 以 最 优 个 体 引导 进 
化 能 加 快 算法 的 收敛 速度 ， 但 不 利于 维持 群体 的 多 样 性 ， 
致 算法 陷入 早熟 收敛 。 多 群 操作 是 一 种 改善 群体 多 样 性 的 有 效 
方法 HH， 子 群 间 信 息 的 交互 能 为 个 体 学 习 提 供 更 多 的 信息 。 多 
群 SGO 对 基本 SGO 算法 的 个 体位 置 更 新 方法 进行 改进 ， 充 分 
利用 子 群 的 最 优 个 体 信息 和 整个 群体 的 最 优 信息 ， 引 入 量子 学 


与 全 局 最 优 解 和 本 群 中 其 他 个 体 信息 相 结 合 ， 增 加 子 群 个 体 学 
习 方 法 的 多 样 性 ;，c) 采 用 简单 的 子 群生 成 方法 ， 既 保证 了 子 群 
个 体 的 充分 进化 ， 又 达到 子 群 间 信息 交换 和 改善 子 群 多 样 性 的 
目的 。 


1 SGO 算 法 


复杂 问题 的 解决 需要 多 方面 因素 的 参与 ， 当 某 个 个 体 不 足 
以 解决 复杂 问题 时 ， 需 要 向 周围 或 拥有 此 项 能 力 的 个 体 学 习 ， 
以 便 提 高 自身 解决 问题 的 能 力 。SGO 算法 是 基于 此 种 思想 而 提 
出 的 一 种 新 型 优化 算法 。SGO 算法 实现 过 程 主要 包含 初始 化 、 
提高 、 获 得 阶段 三 个 步骤 。 描 述 如 下 : 
1.1 种 群 初始 化 

随机 初始 化 种 群 P， 如 式 (1) 所 示 。 

Xij = Kin +r * (1, D) * (Xmax —Xui) (1) 

其 中 : PL2,3,-N, jzL2,3,.—D ; 为 种 群 的 规模 ;，DD 为 变 
量 的 维 数 ;， 7 为 0 到 1 之 间 的 随机 数 ， Xow 与 Xs 分 别 为 变量 
的 上 下 限 ; P 为 进化 种 群 ，X 为 个 体位 置 
12 提高 阶段 


SGO 算法 在 此 阶段 采 | 


的 学 习 策略 如 式 (2) 所 示 。 


Xnew, , = c*Xold, , + r*( gbest, 一 Xold, ,) (2) 


习 方 法 ， 随 机 选择 部 分 个 体 进 行 更 新 ， 以 提高 算法 的 收敛 精度 
和 增加 群体 多 样 性 ， 减 小 算法 陷入 局 部 最 优 的 可 能 。 

2.1 提高 阶段 的 改进 
式 (2) 可 知 ， 基 本 SGO 算法 的 提高 阶段 ， 个 体 向 群体 
的 最 优 个 体 学 习 ， 目 的 是 提高 算法 的 收敛 速度 ， 而 快速 的 收敛 
易 导 致 群体 多 样 性 丢失 ， 算 法 易于 陷入 局 部 最 优 。 此 外 ， 从 知 
识 的 学 习 难 易 程 度 看 ， 个 体 更 容易 从 其 较 近 的 邻 域 个 体 获得 知 
识 。 基 于 这 样 的 思想 ， 在 MPSGO 算法 中 的 提高 阶段 ， 以 子 群 
最 优 个 体 的 信息 对 个 体 进行 引导 ， 利 用 就 近 获取 知识 方法 使 个 
体 在 其 局 部 区 域 增长 知识 ， 增 加 个 体 学 习 知 识 的 多 样 性 。 具 体 
BRE TR: 首先 将 所 有 个 体 随 机 分 成 个 体 数 量 相等 的 若干 个 子 
群 ， 再 求 取 每 个 子 群 的 最 好 个 体 ， 按 照 式 (5) 对 各 子 群 个 体位 
进行 更 新 。 


Xnew, , — c*Xold, ;+ r*(agbest(a, j)- Xold, ,) (5) 


其 中 : c 的 取 值 与 基本 SGO 算法 相同 ，7 为 0~1 间 的 随机 数 ; 
a 表示 第 a 个子 群 ，Xnew;， 与 Xold,, 分别 为 第 i 个 个 体 更 新 
前 后 的 第 j 维基 因 ; agbest(a. j) 为 第 a 个 子 群 中 最 好 个 体 的 第 
j 维基 因 。 
2.29 获得 阶段 的 改进 

X (3) (4) 可 知 ， 在 SGO 算法 的 获得 阶段 ， 个 体 向 群 
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体 中 的 其 他 个 体 及 群体 中 的 最 好 个 体 学 习 以 获 得 新 知识 ， 但 
未 考虑 个 体 学 习 能 力 的 局 限 性 和 学 习 方 法 的 多 样 性 ， 不 利于 群 
体 多 样 性 的 维持 。 基 于 这 样 的 分 析 ， 在 MPSGO 算法 的 获得 阶 
段 ， 被 选择 用 于 学 习 的 随机 个 体 限制 在 子 群 内 部 ， 个 体 除了 向 
子 群 内 其 他 个 体 和 种 群 最 优 个 体 学 习 外 ， 还 增加 了 个 体 向 其 他 
随机 选择 子 群 的 最 好 个 体 学 习 , 以 提高 个 体 学 习 方法 的 多 样 性 。 
改进 的 获得 阶段 个 体 更 新 方程 如 式 (6)(7) Bran. 

对 某 子 群 中 的 个 体 i, 在 其 所 在 的 子 群 中 随机 选择 个 体 k 作 


2.4 子 群 信息 交换 

子 群 间 信息 交换 是 保持 各 子 群 多 样 性 的 一 个 重要 方法 。1 
于 子 群 规模 往往 较 小 ， 当 迭代 到 一 定 代数 时 ， 子 群 中 的 个 体 易 
于 趋同 ， 导 致 子 群 多 样 性 减 小 ， 不 利于 算法 的 收敛 。 此 外 ， 为 
给 子 群 中 的 个 体 充 分 的 进化 时 间 , 各 子 群 必须 进化 一 定 的 代数 。 
为 不 增加 算法 的 复杂 性 ， 在 MPSGO 算法 中 ， 采 用 一 种 简单 的 
分 群 方法 ， 当 群体 进化 一 定 代数 后 ， 所 有 子 群 的 个 体 随 机 重新 
组 合 ， 形 成 新 的 子 群 。 有 具体 操作 方法 为 : 如 果 进 化 代数 能 被 给 


为 其 学 习 的 第 一 部 分 ， 在 所 有 子 群 中 ， 随 机 选择 一 个 子 群 的 最 
优 解 作为 其 学 习 的 第 三 部 分 ， 种 群 的 当前 最 优 解 作为 其 学 习 的 
第 三 部 分 ， 新 的 学 习 机 制 如 下 。 

If X, isbetterthan X, 


Xnew, ;= X, tn «(X - X, )+n *(gbest; -X,) 十 


Hj 
(6 
A *(agbest(aa, j) — X, ;) ) 


else 


Xnew, ;= X utn *(x, , =X; )* n * (gbest, =X; )+ d 
n *(agbest(aa. j) -X,) 


其 中 : aa 为 随机 选择 子 群 ， Xij 与 Xi,; 分 别 为 个 体 i 和 个 体 k 


定 的 子 群 进化 周期 整除 ， 则 所 有 群体 的 个 体重 新 混合 ， 随 机 形 
成 新 的 子 群 (If mod(gen,7)=0 then regroup the population). gen 
为 当前 的 优化 代数 。7 为 设 定 的 整数 值 ，7 越 大 越 有 利于 子 群 
充分 进化 ,，7 越 小 越 有 利于 维持 子 群 多 样 性 。 
2.5 算法 收敛 性 分 析 

最 小 值 优化 问题 如 下 : 

min{ f(X)| VX €Q,0« f (X) «4] 

其 中 Q 是 可 行 搜索 空间 且 是 非 空 有 限 集 。 

假设 1 搜索 空间 是 范围 为 的 DD 维 空间 , 其 中 有 N 个 局 
部 最 优 解 ， 相 对 应 的 区 域 范 围 记 是 5,5,…,5, , 记 KEJI p 
域 S* 的 最 好 解 。 


Ix 


的 第 j 维基 因 ; abest, 为 所 有 子 群 中 最 优 个 体 的 第 j 维基 因 ; 
agbest(aa, j) 为 子 群 aa 中 最 优 个 体 的 第 j 维基 因 ; Xnew,, 为 
Xij 为 个 体 i 更 新 前 后 第 j 维基 因 的 位 置 。 

X (6) CD 表明 , 在 MPSGO 算法 的 获得 阶段 , 个 体 不 仅 
从 其 所 在 子 群 中 获得 知识 ， 同 时 也 从 全 体 子 群 的 最 优 解 和 其 他 
子 群 的 最 优 解 中 获得 知识 ， 既 兼顾 了 个 体 学 习 范围 的 局 限 性 ， 
又 考虑 到 了 个 体 学 习 方 法 的 多 样 性 。 
2.3 ”量子 学 习 策 略 

Grover 在 1996 年 提出 了 一 种 高 效 的 量子 算法 00 ， 即 量子 
搜索 算法 。 随 后 ， 研 究 者 把 量子 计算 与 机 器 学 习 结 合 起 来 ， 形 
成 了 量子 学 习 02 巩 。 量 子 算法 特点 就 是 利用 量子 态 的 干涉 作用 ， 
使 所 需 的 结果 增强 ， 不 需要 的 结果 减弱 ， 这 样 在 检测 时 所 需 的 
结果 就 会 以 较 高 的 概率 出 现 03。MPSGO 算法 中 ， 在 子 群 更 新 
过 程 中 引入 量子 学 习 ， 目 的 是 提高 较 优 解 的 出 现 概率 。 本 方法 
中 ， 以 子 群 最 优 解 (agbest) 及 个 体 与 群体 平均 值 之 差 
C X, - aMean(a) ) 作为 个 体 学 习 的 向 导 。 为 减少 适应 度 函数 的 
评价 次 数 ， 本 文 仅 对 子 群 中 随机 选择 的 部 分 个 体 进行 操作 。 方 


TT 


法 如 下 。 
Xnew, , = temp 十 ((—1) * ceil(0.2 + rand (1, j))) * 
b*(-log(u)) i 
temp =( p.* X, , - m.*agbest(a, j))/ (p +m) (9) 


b — rand (1) * (aMeanX (a, j)- X, ;) (10) 


其 中 : p.m,u 3 0-1 间 的 随机 数 ; 
个 体 的 位 置 平均 解 。 


aMean(.) HTE a 中 所 有 


假设 2 存在 入 个 子 群 , 搜索 区 域 为 4, 2……9, ,每 个 子 群 
可 以 同时 搜索 相对 应 的 区 域 。 

定理 1 如 果 前 两 个 假设 满足 ， 则 多 子 群 算法 比 全 群体 算 
法 具有 更 高 找到 全 局 最 优 解 的 概率 。 


证 明 对 于 全 群 算法 ， 令 p{Xe5*} 为 找到 最 优 解 的 概率 ， 
p{X e S*} 定义 如 下 : 


N N S 
px esy-[Ioix es) =] fẹ a1) 
n=1 


对 多 子 群 算法 ， 令 PiX e5*) 为 找到 最 优 解 的 概率 ， 
Pw{XeS*} 定义 如 下 : 


N 


N 
py{X es*}=[ [p{X eS.) -JJS 


a2) 
n-l n-l n 

N S, 
ji PAX e s*) Dp. NO 13 
^ pes) [T S, 529, zi 

mO 

KAQ «O (-L2, .., N), Bil 

PAX € S*} > p(X e $*] (14) 


上 述 分 析 可 得 ， 多 子 群 MPSGO 算法 找到 最 优 解 的 概率 
优 于 单一 群体 的 SGO 算法 。 
对 于 MPSGO 算法 而 言 ， 群 体 中 个 体 学 习 可 以 认为 是 变异 
操作 算 子 (MO), 子 代 学 习 者 与 父 代 学 习 者 的 比较 可 以 看 做 是 选 
TEBRTE TET (SOUS, 在 每 次 迭代 中 , MPSGO 算法 优化 过 程 依次 
MO 与 SO 实现 ， 故 利用 随机 泛 函 理论 可 以 将 一 次 迭代 抽象 
为 由 MO 和 SO 合成 的 随机 映射 。 
定义 1 变异 算 子 MO) 是 一 种 以 一 定 概率 对 个 体 向 量 的 
各 维度 分 量 进 行 重 组 变换 的 过 程 。 它 是 解 空间 上 的 一 种 随机 映 
JT: MxQo Q" 。 在 MPSGO 中 ， 这 变异 操作 定义 为 
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pio|T, (o. X,) 7 newX,} 
-p(newX, = c*Xold, (15) 
+r*(agbest(a, j) 一 Xold,)) 
其 中 : @ 表 示 基 本 事件 ， 其 他 参数 含义 与 式 (5) 相同。 
定义 2 选择 算 子 (S0) Pd ERES. 
来 选择 子 代 学 习 者 的 过 程 ， 它 是 解 空间 oV 与 解 空 间 Q 之 间 的 


HKR, ARENT, o O >Q. 
对 于 MPSGO 算法 ， 其 选择 算 子 定义 为 
T, (X, newX,) 
-Z^AZeiX,newX,) ^ f (Z) (16) 
-min( f ( X,), f (newX,)] 

由 于 变异 算 子 和 选择 算 子 都 是 随机 生成 的 ， 所 以 产生 的 学 
习 者 也 是 随机 的 ， 使 算法 具有 跳出 局 部 最 优 的 能 力 ， 不 易 陷 入 
局 部 最 优 。 在 MPSGO 的 一 次 迭代 中 ， 相 当 于 将 映射 了 作用 了 
当前 群体 P，7 表 达 式 如 下 所 示 : 

T= T, :MxQ>Q (17) 

4 XO 为 第 ! 次 迭代 时 群体 中 的 第 i 个 个 体 ，X,(t+D 为 第 
1+1 次 友 代 时 群体 中 的 第 :个 个 体 ， 那 么 


| 
| 
H 


X,t*D 2T(o, X, (0) - T, (T,, (c. x,())) (18) 


f£ T 作用 下 ，MPSGO 每 次 迭代 所 产生 的 子 代 个 体 并 不 差 
于 其 相应 的 父 代 个 体 ， 因 此 ， 各 代 群 体 中 最 优 学 习 者 对 应 的 适 
应 值 形 成 一 个 单调 非 增 序 列 ， 即 


 F(X4,(r-1))2 


g| 38 173] d 是 距离 函数 4(X.Y) H (X) - f(,)1 ， 
YX Y, EQ ZLTA QxQ 8| R IB, MWA, (Ld) Ei 
完整 的 赋 范 空间 

定义 3 令 Q 和 R 是 两 个 赋 范 空 
如 存在 一 个 随机 变量 Ko <1 使 得 

pio|d(T(o. X,).T(o. X,..)) 
x K(o)a(r(o.X, ,),T(v. X,)))21 
称 映射 7:o >R 是 一 个 随机 压缩 映射 。 

定理 2 MPSGO 的 一 次 迭代 所 形成 的 随机 映射 7 是 一 个 
随机 压缩 算 子 。 

证 明 据 MO 算 子 与 SO 算 子 的 定义 及 式 (18) 可 知 ， 
MPSGO 在 每 次 欠 代 过 程 中 ， 子 代 最 好 个 体 可 能 优 于 父 代 最 好 
个 体 。 因 而 存在 一 个 随机 变量 Ko <1 ， 使 得 

d(T(o.X, ).T(o.x )) 24(X,. X.) M f(x) -f(x,.)ls 
K(o)| f(x.) - fGc)Ie K(o)a(x, ,.X,) 


F(X pa (D) (Xna lt +) 19) 


间 。 定 义 X,=T(w,X,,)， 


(20) 


Q1) 


. 2 -(o|d(T(o,X,).T(o.X,.)) 
所 以 i uus Ed (22) 


因此 ， p(@)=1 (23) 
在 MPSGO 中 由 MO 和 SO 合成 的 映射 7 是 一 个 随机 压缩 


映射 。 


引 理 207 设 Q 是 一 个 巴 拿 赫 空间 。 若 7T:o~>@ 是 一 个 压 
缩 映 射 ， 则 7 在 Q 内 有 唯一 不 动 点 。 

定理 3 Wi T:MxoooJéÉ MPSGO 中 的 压缩 映射 ， 那 么 ， 
Eo 中 7 有 唯一 不 动 点 ， 也 即 ，MPSGO 是 渐 近 收敛 。 
2.6 SGO 与 MPSGO 的 多 样 性 分 析 
本 文选 取 不 同类 型 的 函数 F2、F4、F6 8I F8 来 检验 算法 的 
多 样 性 。 测 量 多 样 性 的 方法 很 多 ， 如 基于 所 有 个 体 平均 汉 明 距 
离 的 方法 MAMER ELSE FoU IR PS 7r V UC 73 1 
08。 为 了 简单 直观 分 析 SGO 和 MPSGO 算法 的 多 样 性 ， 本 文 
采用 基于 平均 距离 来 衡量 群体 多 样 性 。 为 了 公平 比较 SGO 和 
MPSGO 算法 的 多 样 性 ， 两 个 个 体 之 间 的 距离 被 归 一 化 到 [0,1] 
内 。 群 体 多 样 性 定义 如 下 : 

定义 4 ”群体 中 个 体 p; 5s pi/ 之 间 的 距离 为 


H(p,,p,)= 3 -pl (24) 


定义 5 群体 中 个 体 p; 与 Pp, 之 间 的 归 一 化 距离 为 


H` (psp,)=H(p,,p,)/ max{H (p,,, Pa)» 


Q3) 
m=1,2,...N -Ln=m+l,m+2,...,N} 
定义 6 每 次 群体 迭代 中 的 平均 距离 为 
ary - 2o 2. Gn) Q9 


N*(N -1) 

在 本 次 测试 中 ，SGO 和 MPSGO 算法 的 参数 设置 如 下 : 两 
算法 的 种 群 大 小 均 为 50， 以 算法 最 大 迭代 次 数 作为 终止 条 件 ， 
算法 代数 为 4 000， 函 数 维 数 为 50， 两 种 算法 的 多 样 性 变化 如 
图 1 所 示 。 


sco 
MPSGO. 


i [r3 
8 — 180 — 240 — 320 — 4X0 0 wo — (60 — 200 3200 400 


78 


图 1 SGO 和 MPSGO 的 多 样 性 变化 


1 结果 显示 ， 对 选取 的 四 个 函数 ，MPSGO 的 多 样 性 优 
于 基本 SGO. F2 函数 图 像 表 明 在 初始 化 阶段 ， MPSGO 的 多 样 
性 差 于 基本 SGO, 但 在 中 后 期 ,MPSGO 的 多 样 性 优 于 基本 SGO; 
F8 函数 图 像 显 示 ， 在 迭代 1200-1500 代 时 ，MPSGO 的 多 样 性 
差 于 基本 SGO， 但 在 1 500 代 后 ，MPSGO 的 多 样 性 明显 优 于 
基本 SGO. F4 和 F6 函数 图 像 表 明 ，MPSGO 的 多 样 性 始终 优 
于 基本 SGO 算法 。 
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2.7 ”改进 算法 流程 


以 上 分 析 ，MPSGO 算法 流程 如 下 


初始 化 实验 参数 : 
N,D,maxFEs,maxgen,class,xmin,xmax, 终 


止 准则 


i 


按照 式 (1) 初 始 化 种 群 并 把 种 群 分 组 。 


v 


对 初始 化 种 群 进行 评价 ， 确 定 全 局 最 优 解 


L 


图 2 所 示 : 


提高 阶段 : 按 式 (5) 更 新 种 群 ， 并 对 个 体 
进行 评价 ， 


更 新 当前 解 。 


N Y 
新 解 优 于 当前 解 ? 


i ———— 


确定 各 子 群 中 个 体 及 全 局 最 优 解 


了 


获得 阶段 : 按 式 6,7) 更 新 种 群 ， 并 对 个 体 
进行 评价 ， 


更 新 当前 解 。 


N Y 
ED 新 解 优 于 当前 解 ? 接受 


确定 各 子 群 中 个 体 及 全 局 最 优 解 


Y 


量子 学 习 阶 段 : BOX 6) 进 行 种 群 学 习 ， 
并 对 个 体 进行 评价 ， 更 新 当 解 。 


弃 除 


3 ”仿真 实验 


为 验证 
真实 验 ， 结 果 与 
算法 PH、 
31 参数 设置 

本 文 所 用 的 


F3、F4 为 多 峰 函 


为 50， 


其 他 四 种 算法 参 


下 所 示 : 


PSOFDR 算 法 Pu: 


图 2 MPSGO 算法 流程 


MPSGO 算法 的 性 能 ， 本 文 对 8 个 


NSTLBO 算法 [91、 


SGO 算法 进行 比较 。 


8 ^B 


ETLBO YAPA, 


函数 如 表 1 所 示 。F1、F2 是 单 


基准 函数 进行 仿 
PSOFDR 


Æ K Z, 


数 ，F5~F8 是 旋转 函数 。 所 


D 


算法 的 群体 规模 
以 函数 值 作为 算法 适应 度 值 ， 最 大 函数 评价 次 数 
(maxFES=5000*D) 作 为 算法 的 停止 准则 ， 
群 数量 为 10，c=0.2。 


变量 的 维 数 。 子 


数 均 来 自 相 对 应 的 参考 文 


献 ， 相 关 参 数 如 


93-2; 
NSTLBO $1349: N (neighborhood Size) =3; 
ETLBO 算法 Po]: 


elite size-2; 


SGO THESE c=0.2 
表 1 测试 函数 
问题 函数 变量 范围 最 优 值 
F1 Zakharov [710,10] 0 
F2 Rosenbrock [-2.048,2.048] 0 
F3 Ackley [-32.768,32.768] 0 
F4 Schwefel [-500,500] 0 
F5 Rotated Zakharov [710,10] 0 
Rotated 
F6 [-2.048,2.048] 0 
Rosenbrock 
Rotated 
F7 [-32.768,32.768] 0 
Ackley 
Rotated 
F8 [-500,500] 0 
Schwefel 
3.2 ”实验 结果 及 分 析 
为 验证 算法 解决 不 同 维 优化 问题 的 能 力 , 分别 对 30. 50 维 


函数 进行 实验 。 每 种 算法 独立 运行 30 次 ， TIONS ADS 
的 最 优 解 (best)、 平 均值 (mean)、 标 准 差 (std) 进行 统计 分 
析 ， 实 验 结果 如 表 2、3 所 示 。 表 2、3 中 粗 体 为 最 好 结果 。 
从 表 2 的 结果 可 以 看 出 ,30 维 函 数 条 件 下 ， 对 于 函数 Fl1、 
F3, F5, F7, MPSGO 与 SGO 的 性 能 相当 ， 且 三 项 性 能 都 优 于 
其 他 三 种 算法 ; 对 函数 F2、F4、F8，MPSGO 的 性 能 优 于 其 他 
算法 ， 对 函数 F2、F6，NSTLBO 的 标准 差 优 于 其 他 四 种 算法 ; 
五 种 算法 在 函数 F3 的 标准 差 相同 。 统计 结果 表明 ，MPSGO 的 
平均 性 能 优 于 其 他 算法 。 
从 表 3 的 结果 可 以 看 出 , 在 50 维 函 数 实 验 中 , 对 函数 Fl、 
F3、F5 fll F7, MPSGO 与 SGO 的 性 能 相当 ， 且 三 项 性 能 都 优 
于 其 他 三 种 算法 ; 对 函数 F4 中 ,MPSGO 三 种 性 能 优 于 其 他 算 
法 ; fEF6. F8 实验 结果 中 ，MPSGO 的 平均 值 优 于 其 他 算法 ; 
TE F8 函数 实验 中 ,NSTLBO 算法 的 标准 差 优 于 其 他 四 种 算法 ; 
在 对 F3、F7 函数 实验 中 ，NSTLBO、MPSGO、SGO 三 种 算法 
标准 差 相同 ; 在 函数 F2 实验 结果 中 , PSOFDR 的 最 优 解 、 平 均 
值 优 于 其 他 算法 。 2. 3 统计 分 析 可 看 出 ，MPSGO 的 平均 
性 能 优 于 其 他 算法 ， 尤 其 在 高 维度 问题 下 ， 其 优势 较 明显 。 
为 了 更 直观 地 显示 各 算法 平均 适应 度 的 变化 过 程 ， 本 文 给 
出 了 不 同 算法 对 50 维 函 数 实验 的 优化 过 程 中 平均 适应 度 变化 ， 
如 图 3 所 示 。 图 3 显示 , 对 函数 F3、F4、F6、F7、F8, MPSGO 
算法 的 收敛 速度 快 于 其 他 算法 ; MPSGO 在 对 函数 F2. 实验 前 期 
收敛 速度 快 于 其 他 算法 ， 后 期 收敛 的 精度 差 于 PSOFDR， 但 比 
SGO 算法 效果 好 。 虽 然 在 个 别 函数 优化 实验 中 ，MPSGO 的 收 
敛 速度 不 是 最 快 ， 但 对 于 复杂 函数 ， 特 别 是 旋转 函数 ，MPSGO 


FUP 


comin = 0.4, 0 max —0.9,91 21,92 —1, 


表现 出 了 良好 的 性 能 


录用 稿 刘 亚 军 ， 等 : 


为 了 进一步 证 明 MPSGO 算法 与 其 他 算法 的 差异 ， 采 用 上 
- test 方法 对 各 种 算法 的 结果 进行 测试 。 本 文采 用 双边 检验 ， 显 

性 能 指标 PSOFDR NSTLBO ETLBO SGO MPSGO n : E . 
KFN 0.05. Æ 30. 50 维 条 件 下 ，MPSGO 算法 对 于 其 他 算 

best 1.99E-09 3.30E-73 1.93E-36 0.00E+00 0.00E+00 
y l z= Q? TRES eZ 
Fl mean 1.57E-08 1.03E-70 2.00E-34 0.00E-00 0.00E+00 法 的 + 值 和 p 值 如 表 4、5 所 示 。 表 4、5 ps B, S, W3 
std 1.86E-08 1.65E-70 4.41E-34 0.00E+00 0.00E+00 别 表 示 MPSGO 算法 性 能 优 于 、 等 于 、 劣 于 其 他 算法 的 数量 。 
best 1.19E+01 2.40E+01 1.55E401 2.79E+01 1.29E+01 MPSGO 与 其 他 算法 之 间 的 较 好 结果 以 粗 体 显示 。 由 表 4、5 果 
F2 mean 1.41E+01 245E«01 1.72E+01 2.83E+01 1.34E+01 可 以 看 出 ， 在 采用 双边 检验 ， 显 著 水 平 0.05 情况 下 ，MPSGO 
std 2.41E+00 5.07E-01 1.19E+00 2.49E-01 4.26E-01 x 


数 算法 ， 解 可 以 被 接受 。 
best 7.]1E-15 3.55E-15 3.55E-15 0.00E+00 0.00E+00 多 数 算 ; 解 可 以 被 接受 


表 2 30 维 函 数 测试 结果 


出 


F3 mean 7.11E-15 3.55E-15 3.55E-15 0.00E+00 0.00E+00 a n 008 30748 s. 08 
std 0.00E+00 0.00E-00 0.00E+00 0.00E-00 0.00E-00 denm un, y i 
best 2.82E+03 3.07E+03 3.59E+03 4.97E+03 1.78E+03 E Mmm I 
F4 mean 3.36E+03 443E*03 4.02E*03 5.61E+03 2.40E403 lw 到 |. 
std 4.71E+02 8.11E+02 4.83E+02 6.06E+02 4.35E+02 Pa M H 
best 6.12E-10 2.60E-75 1.84E-38 0.00E+00 0.00E+00 CR | E m 
F5 mean $5.55E-09 5.94E-71 7.80E-36 0.00E+00 0.00E+00 z | | e ， 
std 7.19E-09 1.31E-70 1.04E-35 0.00E+00 0.00E+00 PES l | "| I owe 
_ best — 245E*01 3.13E+01 1.62E+01 1.29E+02 1.17E+01 mi a T Taf RE Dn 
F6 mean  6.83E*01 5.69E+01 3.63E+01 4.16E402 3.70E+01 ep NEED NU | 
std 5.46E+01 2.28E+01 2.44E+01 2.29E+02 2.58E+01 | Di — nro 4 E 
Enc | 
best 7.11E-15 3.55E-15 3.55E-15 0.00E+00 0.00E+00 iso | woo 


F z 
-6 box =H 


F7 mean 1.86E-01 3.55E-15 3.55E-15 0.00E+00 0.00E+00 
Ta std 4.16E-01 0.00E+00 0.00E+00 0.00E+00 0.00E+00 


-8 1 d 


| 
log, (Mean Fitness) 


10 

best 1.84E+03 4.15E+03 3.16E+03 5.27E*03 2.01E+03 24 m S 
F8 mean 2.68E+03 4.92E*03 3.88E+03 6.47E+03 2.48E+03 | al - - | 
std 7.38E+02 7.98E+02 4.69E+02 1.09E+03 4.33E+02 B 05 1 15 2 25 My 08 0 1 2 — ds 
FE 5 FE 5 
. " | 10 l s 10 
一 一 PSOFDR E 
[m — —NSTLBO,. | ， NSTLBO 
表 3 50 维 函 数 测试 结果 tee M guo | 
50 - —-— MPSGO [| | : nl 


函数 性 能 指标 PSOFDR NSTLBO ETLBO SGO MPSGO 
best — 5.28E-04 1.08E-70 2.33E-31 0.00E+00 0.00E+00 
Fl mean  2.64E-02 4.27E-68 2.96E-30 0.00E+00 0.00E-00 
std — 3.29E-02 6.14E-68 3.86E-30 0.00E+00 0.00E+00 
best — 2.67E*01 4.33E+01 3.63E+01 4.80E+01 2.75E+01 | | 
F2 mean 2.90E+01 4.38E+01 3.73E+01 4.85E+01 2.95E+01 ow A m E 
std — 1.39E«00 3.33E-01 6.11E-01 3.46E-01 1.29E+00 
best 2.13E-14 3.55E-15 3.55E-15 0.00E+00 0.00E+00 
F3 mean 3.29E-01 3.55E-15 3.55E-15 0.00E+00 0.00E-00 
std — 7.36E-01 0.00E-00 0.00E-00 0.00E+00 0.00E+00 
best — 446E*03 7.29E+03 6.40E+03 9.14E+03 4.05E+03 
F4 mean 5.89E+03 8.01E+03 7.98E+03 1.01E+04 4.67E+03 


logioCMean Fitness) 


45r ————À - 
-  PSOFDR = 一 PSOFDR 
-— NSTLBO+ NSTLBO 
ETLBO ETLBO 
SGO SGO | 
MPSGO MPSGO 


log, (Mean Fitness) 


std — 1.03E+03 7.00E-02 1.07E-03 7.00E-02 5.94E-02 du. [ | a X mpm 

best — 2.34E-03 3.75E-74 9.18E-34 0.00E+00 0.00E+00 aru Te 
F5 mean 1.40E-02 4.82E-68 6.95E-31 0.00E+00 0.00E+00 = xm ul x 

std — L33E-02 9.86E-68 1.46E-30 0.00E+00 0.00E-00 图 3 8 个 函数 平均 适应 度 值 变化 图 

best 5.12E+01 3.83E+01 4.93E+01 6.79E402 2.23E+01 
F6 mean 7.26E+01 1.07E+02 8.75E+01 1.66E+03 5.68E+01 4 ”结束 语 

std 1.94E+01 4.83E+01 2.14E+01 1.24E+03 3.06E+01 

best 1.03E+00 3.55E-15 3.55E-15 0.00E+00 0.00E+00 本 文 针 对 SGO 算法 收敛 精度 低 及 群体 多 样 性 易 丢 失 问 题 ， 


F7 mean 1.32E*00 3.55E-15 4.26E-15 0.00E-00 0.00E-00 提 
std 2.07E-01 0.00E+00 1.59E-15 0.00E+00 0.00E+00 群 
能 
有 


出 了 基于 多 子 群 的 社会 群体 算法 (MPSGO )。 通 过 引入 多 子 
量子 学 习 及 子 群 信息 交换 方法 , 提高 了 SGO 算法 的 整体 性 
。 通 过 与 SGO 及 其 他 三 种 算法 比较 实验 , 验证 了 改进 方法 的 


best 5.68E*03 8.96E+03 8.38E+03 1.16E+04 3.94E+03 
F8 mean 6.84E+03 9.81E+03 9.61E+03 1.26E+04 5.26E+03 
std 1.20E+03 8.15E*02 9.43E+02 8.88E+02 8.25E+02 


录用 入 
表 4 各 算法 对 8 个 基准 函数 为 30 ZEIT] t-test 结 
函数 指标 PSOFDR NSTLBO ETLBO SGO 
T 4.6258 3.4100 2.4821 0.0000 
P 0.0000 0.0012 0.0156 0.0000 
d T 1.6905 92.2834 16.4651 165.6072 
P 0.0963 0.0000 0.0000 0.0000 
E T 21.5191 13.7595 13.7595 0.0000 
P 0.0000 0.0000 0.0000 0.0000 
" T 8.1838 12.0438 13.6443 23.5496 
P 0.0000 0.0000 0.0000 0.0000 
" T 4.2326 2.4899 4.1208 0.0000 
P 0.0001 0.0157 0.0001 0.0000 
" T 2.8442 3.1724 -0.1087 9.0180 
P 0.0061 0.0024 0.9141 0.0000 
" T 2.4495 13.7595 13.7595 0.0000 
P 0.0174 0.0000 0.0000 0.0000 
M T 1.2860 14.7259 12.0255 18.5730 
P 0.2053 0.0000 0.0000 0.0000 
B 6 8 7 4 
S 2 0 1 4 
W 0 0 0 0 
表 5 各 算法 对 8 个 基准 函数 为 50 维 的 t-test 结 
函数 “指标 PSOFDR NSTLBO ETLBO SGO 
T 5.6831 4.9159 5.4326 0.0000 
s P 0.0000 0.0000 0.0000 0.0000 
" T -1.9850 75.6295 38.6010 100.3924 
P 0.0519 0.0000 0.0000 0.0000 
T 3.1623 13.7595 13.7595 0.0000 
"n P 0.0025 0.0000 0.0000 0.0000 
m T 7.2875 25.7462 19.1339 42.1701 
P 0.0000 0.0000 0.0000 0.0000 
T 7.4309 3.4557 3.3732 0.0000 
i P 0.0000 0.0010 0.0013 0.0000 
T 3.0693 6.2383 5.8008 9.1123 
i P 0.0033 0.0000 0.0000 0.0000 
T 45.0210 13.7595 18.9741 0.0000 
i P 0.0000 0.0000 0.0000 0.0000 
4 T 7.6633 27.7548 24.5591 42.6334 
P 0.0000 0.0000 0.0000 0.0000 
B 7 8 8 4 
S 1 0 0 4 
W 0 0 0 0 
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